Nonuniform Self-Organized Dynamical States in Superconductors with Periodic 

Pinning 



O 
O 

(N 



(N 
O 

o 



^3 

o 
o 



in 
o 

O 

-i— > 



i 

c 

o 
o 



X 



Vyacheslav R. Misko 1,2,3 , Sergey Savel'ev 1,4 , Alexander L. Rakhmanov 1,5 , and Franco Nori 1,2 
Frontier Research System, The Institute of Physical and Chemical Research (RIKEN), Wako-shi, Saitama, 351-0198, Japan 
2 MCTP, Department of Physics, University of Michigan, Ann Arbor, MI 48109-1040, USA 
3 Department of Physics, University of Antwerpen (CGB), B-2020 Antwerpen, Belgium 
4 Department of Physics, Loughborough University, Loughborough LE11 3TU, United Kingdom and 
5 Institute for Theoretical and Applied Electrodynamics Russian Academy of Sciences, 125412 Moscow, Russia 

(Dated: February 6, 2008) 

We consider magnetic flux moving in superconductors with periodic pinning arrays. We show 
that sample heating by moving vortices produces negative differential resistivity (NDR) of both N 
and S type (i.e., N- and S-shaped) in the voltage-current characteristic (VI curve). The uniform 
flux flow state is unstable in the NDR region of the VI curve. Domain structures appear during the 
NDR part of the VI curve of an N type, while a filamentary instability is observed for the NDR of 
an S type. The simultaneous existence of the NDR of both types gives rise to the appearance of 
striking self-organized (both stationary and non-stationary) two-dimensional dynamical structures. 

PACS numbers: 74.25.Qt 



Semiconductor devices Q exhibiting Negative Differ- 
ential Resistivity (NDR) and Conductivity (NDC) have 
played a very important role in science and technol- 
ogy. These useful devices include yj: Gunn effect 
diodes, pnpn-j unctions, etc. Here we study supercon- 
ducting analogs of these semiconductor devices. Ta- 
ble 1 briefly compares NDR in superconductors, semi- 
conductors, plasmas, and manganites. Conceptually, the 
non-uniform self-organized structures (e.g., filaments and 
overheated domains with higher or lower electric fields) 
are related for superconductors, plasmas, and semicon- 
ductors. However, the physical mechanism giving rise to 
the instability of the homogeneous state can be different 
in each case. 

The magnetic flux behavior in superconductors with 
artificial pinning sites has attracted considerable atten- 
tion due to the possibility of constructing samples with 
desired properties 0, 0, 0, 0, 0] . Among such sys- 
tems, samples with Periodic Arrays of Pinning Sites 
(PAPS) are studied intensely because advanced fabri- 
cation techniques allow to design well-defined periodic 
structures with controlled microscopic pinning parame- 
ters. For such systems, Ref.ll4lh.as revealed the existence 
of several dynamical vortex phases, which are similar to 
the ones shown in the inset of Fig. 1 . At low current den- 
sity j, the vortices are pinned and their average velocity, 
v, is zero (phase I). Here v — N^ 1 ^ ■ x, is the average 
vortex velocity in the x direction. At higher j, intersti- 
tial vortices start to move and the flux velocity increases 
with the drive (phase II). Then, a sharp jump in v(j) 
occurs since a fraction of the vortices depins and a very 
disordered uniform phase arises (phase III). This phase 
is analogous to the uniform electron motion in semicon- 
ductor devices. After increasing the applied current j (for 
superconductors) or the applied voltage (for semiconduc- 
tors), the vortex (electron) velocity shows a remarkable 
and non- intuitive sudden drop (i.e., the velocity drops 
even though the applied force increases) . Indeed, when 



j exceeds a threshold value, only incommensurate vor- 
tex rows move (phase IV) . Finally, the increasing driving 
force completely overcomes the pinning (phase V). The 
dependence v(j) describes the voltage-current character- 
istic ( VI curve) since the electric field E is related to v 
by E = —v(j)B/c. Here we prove that the flux motion in 
samples with PAPS results in an unusual, for supercon- 
ductors, VI curve with NDR of the so-called S type 0,0]. 
That is, within some interval of voltages there exist three 
different current values corresponding to a single elec- 
tric field. Such types of 5-shaped VI curves play a very 
important role in plasmas and semiconductors, and give 
rise to a filamentary instability when a uniform current 
flow breaks into filaments with lower and higher current 
densities |l|. 

Sufficiently strong disorder in the pinning array gives 
rise to the disap pea rance of the described dynamical 
phase transitions [14j and, consequently, the vanishing of 
the NDR for S'-shaped VI curves. Higher thermal fluc- 
tuations also result in such an effect. Experimentally, 
the current density, j, at which the flux flow regime is 
observed in superconductors, is usually high @,EH, and 
the Joule heat should also affect the picture described 
above. An increase in temperature T results in a de- 
crease of the pinning force. Thus, the current density can 
decay with increasing electric field, and a iV-shaped VI 
curve with a NDR of N type (red dashed line in Fig. 1) is 
commonly observed in superconductors for high current 
density HEft 

The uniform state in samples with NDR 
of N type is also unstable [l| ; and a propagating resistive 
state boundary or the formation of resistive domains can 
destroy the uniform superconducting state 0, 0, E] • For 
certain pinning parameters and cooling conditions, one 
can achieve a situation where the NDR of both N and S 
type simultaneously coexist in the VI curve (Fig. 1). In 
this case, we predict remarkable flux flow instabilities. 

Here, we study the effect of Joule heating and disorder 
on the VI curve of superconductors with PAPS. Based 
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NDR 
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TABLE I: Comparison between non-uniform non-equilibrium states in superconductors, semiconductors, plasmas, and man- 
ganites with VI curves having Negative Differential Resistivity (NDR). Since IV curves in semiconductors map to VI curves 
in superconductors, then NDC (for semiconductors) maps into NDR for superconductors. Here, N(S) type shapes for semi- 
conductors correspond to S(N) type for superconductors @]. The Negative Differential Conductivity (NDC) found in Q is 
analogous to the Gunn effect in semiconductors, where electron-charge modulations lead to steps in j(E) in the NDC regime. 

on analytical and numerical analysis of the VI curve, we 
find the conditions under which the VI curve has a NDR 
region of either S type or N type, or both. We discuss 
the effect of the shape of the NDR on the vortex and 
current flow. We argue that the coexistence of the NDR 
of both, S and N types, gives rise to the macroscopically 
non-uniform self-organized dynamical structures in the 
flux flow regime. 

We numerically integrate the two-dimensional over- 
damped equations of motion 0, 0, for the flux lines 
driven, by the current j, in the x direction over a square 
PAPS with lattice constant a: rjVj = F™+Ff+Ff +F d . 
Here 77 is the flux flow viscosity , Vj is the velocity of 
ith vortex, F™ is the force acting on the ith vortex per 
unit length due to the interaction with other vortices, 
F^ p is the ith vortex-pin interaction, Fd = j<f>o/c is the 
driving force, and 4>q is the flux quantum. The ther- 
mal fluctuation contribution to the force, Fj\ satisfies: 




(F?(t)) t = and (F?(t)F?(t)) t = 2 V k B TS tJ 6(t - t<) 



The vortex-vortex interaction is modelled by F™ = 
(</»g/87r 2 A 3 (T))Ef=i^i(|r 4 -r,|/A(T))^ : K x is the 
modified Bessel function, the summation is performed 



over the positions i\- 



in 



,0/|r< 



of N v vortices in the sample, and 
— Tj I . The temperature depen- 



dence of the penetration depth A is approximated as 
A(T) = A (l - T 2 /^)- 1 / 2 . The Ginzburg-Landau for- 
mula for H C 2(T) is used and we assume that the ra- 
tio kgl = A/£ is temperature independent (here £ is 
the coherence length). The pinning is modelled as N p 

parabolic wells located at positions r). . The pinning 

force per unit length is Ff = (F p (T)/r p ) £f=i |r 4 - 
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r[ p ^|/AoJ f^ ; , where r p is the range of 
the pinning potential, is the Heaviside step function, 
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We estimate the max- 

2 £2 , 



imum pinning force, F p — F p (T), as H 2 t; /r p and, thus 
F p (T)^F p0 (l-T 2 /T 2 ). 



FIG. 1: (Color online) The average vortex velocity v oc E 
vs current j for B/B® = 1.074, r p = 0.2A and F p /F = 2 
for increasing (red open circles) and for decreasing (blue open 
squares) j. State A is unstable and the sample divides into 
filaments, some in states B and some in C. State C is also 
unstable. The corresponding stable states are on the lower 
(point E) and on the upper (point D) VI curve branches. 
Inset: VI curve when no heating effects are taken into account. 

For brevity, the simulations shown here are for 18 x 
12 Ag periodic cells and at magnetic fields near the first 
matching field, = cf> /a 2 , where N p — N v . First, 
we start with a high-temperature vortex liquid. Then, 
the temperature is slowly decreased down to T = 0. 
When cooling down, vortices adjust themselves to min- 
imize their energy, simulating field-cooled experiments. 
Then we increase the driving current and compute the 
average vortex velocity v(j). 

The average power of Joule heating per unit volume 
is jE = jvB/c. We assume that the temperature relax- 
ation length is larger than any local scale and the sample 
thickness. Under such conditions, the local temperature 
increase due to vortex motion can be found from the heat 
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balance equation 0: hoS (T — To) = vjBV/c, where ho 
is the heat transfer coefficient, T the ambient tempera- 
ture, and S and V are the sample surface and volume. 
Further, we shall assume that To <C T c and neglect T). 
We define the dimcnsionlcss driving force as /,/ = j /jo 
and introduce the dimcnsionless parameters V x = v/vq 
and b = B/B^, where vq = c 2 / '4wKQ L a n Xo, and jo — 
c^o/87t 2 Aq. We assume that the normal state conductiv- 
ity a n is temperature-independent. As a result, the tem- 
perature and the driving force are related by: T/T c = 
K t hV x f d b, where K th = j v B$V/ch T c S is the ra- 
tio of the characteristic heat release to heat removal. 
Using the rough estimates A = 2000 A, kql = 100, 
cr„ = 10 16 s" 1 , V/S = 1000 A, T c = 90 K, B^ = 500 G, 
and ho = 1 W/cm 2 K, we find K th = 0.05 - 0.06 and F p0 
is of the order of Fq = 4>ojo/c. In the simulations, we 
used K th = 0.0525. 

The velocity v versus current j (which coincides with 
the sample VI curve in dimensionless variables) is pre- 
sented in Fig. 1 for increasing and for decreasing j (if 
we neglect the heating effect, v(j) has the shape shown 
in the inset of Fig. 1, similar to Ref. 14] ). For low cur- 
rents, j < 3, the effect of heating is negligible; for j > 3.7, 
the behavior of v(j) drastically changes, compared to the 
non-heating case shown in the inset of Fig. 1. In particu- 
lar, an abrupt transition occurs between regimes IV and 
V. The most pronounced feature related to heating is the 
appearance of hysteresis in regions IV and V: the over- 
heated vortex lattice (for decreasing j) keeps moving as a 
whole at lower currents than the "cold" one (for increas- 
ing j). As a result, we obtain a complex N and S type 
VI curve, characterized by two kinds of instabilities. For 
example, if the current density exceeds the value j w 3.7 
(point A), the uniform current flow is unstable and the 
so-called filamentary instability Q] occurs. Consequently, 
the current flow breaks into supercurrent filaments, some 
with lower js (state B) and others with higher jc (state 
C). The state C is, in its turn, unstable. The correspond- 
ing stable states are on the lower (E) and on the upper 
(D) VI curve branches. 

Note that a small amount of pinning disorder influ- 
ences v(j) and can lead to the disappearance of phase III 
|19|. However, the robust hysteresis, related to heating, 
remains. 

For a sample included in an electric circuit (see inset 
in Fig. 2), the circuit equation is LI + RI + IE = U, 
where I = jA, L and R are the inductance and resis- 
tance, A and I are the sample cross-section and length, 
and U is a constant voltage. The circuit and Maxwell 
equations describe the development of small field pertur- 
bations S~E and <5B and current Sj. We seek perturba- 
tions of the form SE, SB, Sj cx exp(^i/to), where v is the 
value to be found and to — L/R. An instability develops 
if Re{v) > 0. In general, we should consider also the ther- 
mal equation, but for filamentary instabilities the tem- 
perature rise is not of importance. To find the instability 
criterion we consider y-dependent perturbations p]. In 
such a geometry, SE has only the x component, while <5B 
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FIG. 2: (Color online) Schematic VI curve of the supercon- 
ductor for two different values of the hysteresis due to over- 
heating. The more pronounced hysteresis loop (shown by the 
green solid line) corresponds to a larger value of the charac- 
teristic heat release. Inset: electrical circuit and the sample 
with filaments. 

has only the z component. Using SB = ctoSE' j 'vw, we 
find from the circuit and Maxwell equations that: 



(v+l^j + p-^E = 0, 

5E" - 0SE' - ^SE = 0, 
to 



(1) 
(2) 



where Sj and SE are the values averaged over the cross- 
section, prime is differentiation over the dimensionless 
coordinate y/w, w is the sample half-width, p c — RA/l, 
t s = (ATTw 2 /c 2 )(dj/dE), p = (Airw/c){dj/dB) are pa- 
rameters, which are either positive or negative depending 
on the part of the VI curve for given background fields E 
and B. All quantities are averaged over a volume which 
includes a large number of vortices. 

The first boundary condition to Eq. (J2J, SB(w) = 
—5B(—w), is obtained assuming that the applied mag- 
netic field B is constant. From Maxwell equations we 
derive SB(w) — SB(—w) = 8nwSj/c. Using these SB's, 
and substituting Sj from Eq. ijTJl. we obtain the second 
boundary condition SE'(w) — ~^vSE j {v + 1), where 
7 = 4ttw 2 1/c 2 AL. Substituting the solution of Eq. J5J, 
SE = C\ exp(piy/w) + C2 exp(p2j//w), to the boundary 
conditions and requiring a zero determinant for the ob- 
tained set of linear equations for constants Ci, we find 
the equation for the eigenvalues v 
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{V + 1)P2 



Pi tanhj>2 
P2 tanhpi 



Pi 



71/ 



(3) 



where Pl , 2 = (3/2 ± ^/(3 2 /A + vt s /t Q . 

For small values of \dj/dB\, the solution of Eq. © 
can be found explicitly (with accuracy up to |/3 2 |): v = 
— 1 — jto/t s . It follows from this relation that the insta- 
bility occurs only at the VI curve branch with NDR, if 
the drop of the voltage is large: \dE/dj\ > p c . The char- 
acteristic size of the arising filamentary structure is of the 
order of Aw oc w/\pi \ = w/y/\j + t s /t \. If \^ft s /t \ > 1, 
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the filament width is small, Aw <C w. Thus, the sample 
with NDR in the VI curve of S type divides into small 
filaments with different current densities (in different dy- 
namic flux flow phases III and IV) at Rl/A <C \dE/dj\ 
and L <C 4irl/c 2 A. The obtained results are valid if 
p~ 1 \dE/dj\ » [(ATTw/c)dj/dB] 2 . The instability occurs 
only if the left hand side of this last inequality is higher 
than unity; the right hand side is much smaller than 1 
for the parameter range used above if to < 1 mm. In 
the stationary inhomogeneous state that arises after de- 
velopment of the instability, the electric field should be 
uniform over the sample. 

A more complex dynamics appears when the VI curve 
has simultaneously both NDR parts of N and S types 
(see Fig. 2). In this case, the filaments with higher cur- 
rent density jc are unstable if the system is far from the 
voltage-bias regime P, Q • The instability of the filament 
with an TV type VI curve results in the switch of the fil- 
ament into the overheated state [16| D, giving rise to a 
non-uniform electric field distribution in the sample and, 
as a result, to non-zero B. Thus, the state that appears 
after the instability develops is non-stationary. We con- 
sider two possible VI curves (type 1 red, type 2 green) 
shown in Fig. 2. In state D the flux lines move fast, which 
is accompanied with the acceleration of the flux flow in 
the lower-current filaments. In the nearly current-biased 
mode |l| , if the VI curve is of type 2, the sample comes to 
the uniform stable state A' with the current density ja- 
However, if the VI curve has the form shown in the (red) 
curve 1, the stable uniform state with the current den- 
sity ja does not exist. In this case the high electric field 
state moves from point D to point F and falls down to 
a lower branch of the VI curve (point F'). In this state, 
the electric field is also non-uniform and the lower and 
higher resistivity states will move towards A. However, 



this state is unstable and the cycle of transitions will be 
repeated (A — > C — > D — > F — > F' — > A). Such a cyclical 
dynamic state could be realized in the form of flowing 
either stripes or resistive domain walls moving along the 
filaments. 

It is important to stress that the described cyclical 5- 
step dynamics (A — > C — > D — > F — » F 1 — » A) obtained 
for the IV curve having NDR of both N and S types can- 
not be realized for IV curves with either only N or only S 
type of NDR. The appearance of the non-stationary os- 
cillatory regime for stationary external conditions is very 
unusual. This cyclic dynamics could be extended to dif- 
ferent physical systems, e.g., plasmas or superconductors 
without artificial pins driven by a current flowing along 
the externally applied magnetic field |2jj. Moreover, the 
predicted dynamical behavior (which can be generalized 
for several other systems, e.g., semiconductors) is poten- 
tially useful for the transformation of a dc input into 
either an ac-current or a voltage output which can be 
controlled by the parameters of the external circuit. In 
a broader sense, this general class of cyclic dynamics (ac 
output from dc input) is also found in other important 
nonlinear systems, like the dc Josephson effect. 

In summary, for superconductors with periodic pinning 
arrays with certain pinning and heat transfer characteris- 
tics, we derive VI curves with NDR of N type, S type, or 
both. Complex dynamics and regimes including domain 
structures and filamentary instabilities appear when the 
VI curve has both S and N types of NDR. We ana- 
lyzed the self-organized non-uniform dynamical regimes 
for these instabilities. 

This work was supported in part by ARDA and NSA 
under AFOSR contract F49620-02-1-0334; and by the US 
NSF grant No. EIA-0130383. 
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